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ABSTRACT 

Origins  of  Convective  Activity  Over  Panama.  (August  1989) 
Christopher  Stephen  Strager,  B.S.,  The  Pennsylvania  State  University 
Chair  of  Advisory  Committee:  Dr.  Steven  W.  Lyons 

Satellite-derived  outgoing  longwave  radiation  (OLR)  data  were  used  to  examine 
convective  variability  over  the  Panama  region.  Time  series  analysis  of  the  area- 
averaged  daily  OLR  data  for  1984  and  1985  revealed  a  persistent  12-d  oscillation  in 
convective  activity  during  each  season.  Composite  analyses  of  OLR  data  for  the  area 
120°W-40°W  and  35°S-35°N  for  the  1984  dry  (1  January-9  May)  and  wet  (10  May- 
4  December)  seasons  showed  this  oscillation  extends  beyond  Panama  and  the  Central 
America j^egion  into  the  eastern  Pacific.^ 

Composite  diagrams  of  u-  and  v-wind  components  and  geopotential  heights  were 
constructed  using  daily  upper-air  soundings  from  Albrook  AFS,  Panama  in  order  to 
infer  a  vertical  structure  to  the  oscillation.  These  composites  showed  well-defined 
trends  in  all  three  variables  associated  with  each  phase  of  the  12-d  oscillation. /The 
dry  season  OLR,  wind,  and  geopotential  height  composites  suggest  that  a  west-to- 
east  moving  trough  in  the  upper-levels  seems  to  be  responsible  for  the  convection 
anomalies  in  the  Panama  region.  Wet  season  OLR,  wind,  and  geopotential  height 
composites  suggest  that  periodic  westward  migrations  of  an  upper-level  anticyclone 
associated  with  the  subequatorial  ridge  could  be  producing  the  west-northwestward 
moving  periodic  convective  fluctuations  which  originate  over  Panama  during  that 
season. 
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A  modulation  of  the  12-d  oscillation  was  found  during  the  1984  wet  season  with 
a  period  of  approximately  70-80  d.  It  appears  to  be  connected  with  tropichl  cyclone 
activity  in  the  eastern  Pacific.  For  periods  when  the  12-d  oscillation  was  weak  (well- 
defined)  over  the  Panama  region,  a  significant  increase  (decrease)  in  the  number  of 
tropical  cyclones,  as  well  as  an  increase  (decrease)  in  the  number  of  systems  that 
developed  into  hurricanes,  was  detected* 
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CHAPTER  I 

INTRODUCTION 

The  accurate  forecasting  of  significant  precipitation  events  in  the  tropics  con¬ 
tinues  to  be  a  problem  for  the  operational  forecaster.  Recent  advances  in  satellite  and 
radar  technology  have  allowed  forecasters  to  produce  reasonable  short-range  forecasts, 
but  are  of  little  use  in  developing  medium-  and  long-range  forecasts.  The  scarcity, 
as  well  as  the  questionable  accuracy  of  conventional  data  throughout  the  tropical 
regions,  adds  to  the  problem. 

These  forecasting  problems  are  compounded  in  Panama  where,  in  the  summer 
wet  season  (May- November),  daily  convective  activity  can  be  greatly  enhanced  by  the 
passage  of  waves  in  the  easterlies.  The  northward  migration  of  the  monsoon  trough 
can  also  produce  significant  increases  in  cloud  cover  and  precipitation  amounts.  In 
the  dry  season,  precipitation  events  along  the  north  coast  and  central  highlands 
of  Panama  can  be  initiated  by  the  southward  movement  of  shear  lines  into  the 
Caribbean.  Early  detection  of  these  events  is  often  the  key  to  forecasting  precipitation 
events  beyond  24  h.  In  addition,  spectral  analyses  of  tropical  wind,  pressure,  and 
outgoing  longwave  radiation  (OLR)  fields  over  the  past  20  y  have  revealed  periodic 
fluctuations  in  these  variables  ranging  from  3  to  60  d.  Monitoring  the  progress  of 
these  oscillations  could  be  a  key  to  improving  medium-  to  long-range  forecasting  in 
the  Panama  region  and  in  the  tropics  in  general. 


The  style  used  is  that  of  the  Monthly  Weather  Review 
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In  this  study,  spectral  analysis  techniques  are  applied  to  time  series  of  OLR 
over  the  Panama  region  to  determine  if  any  significant  periodic  fluctuations  in  the 
OLR  field  are  evident.  Next,  composite  OLR  diagrams  are  constructed  to  determine 
the  spatial  extent  and  strength  of  the  periodicities.  Daily  rawinsonde  analyses  are 
also  used  to  infer  the  vertical  structure  and  synoptic  patterns  associated  with  the 
oscillation. 

Chapter  II  contains  a  review  of  previous  work  concerning  tropical  weather 
systems  and  their  time  scales.  The  objectives  and  general  procedures  of  this  research 
are  outlined  in  Chapter  III.  A  description  of  the  data  sets  can  be  found  in  Chapter 
IV.  Chapter  V  describes  the  methods  used  in  the  time  series  analysis  of  the  OLR  and 
rawinsonde  data  sets.  The  results  of  this  research  are  presented  in  Chapter  VI,  along 
with  a  case  study.  Chapter  VII  shows  an  interesting  relationship  between  tropical 
cyclone  activity  in  the  eastern  Pacific  and  convective  variability  over  Panama.  Finally, 
a  summary  and  discussion  are  contained  in  Chapter  VIII. 
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CHAPTER  II 

PREVIOUS  WORK 

Over  the  past  50  y,  many  studies  have  attempted  to  examine  convective  variabil¬ 
ity  in  the  tropics.  These  studies  revealed  oscillations  in  atmospheric  variables  ranging 
from  three  days  to  greater  than  60  d.  The  first  synoptic  feature  detected  which  pro¬ 
duced  significant  convective  variability  over  Panama  were  waves  in  the  easterlies. 
Dunn  (1940)  was  one  of  the  first  to  recognize  these  disturbances  in  the  Caribbean 
region  as  he  tracked  isallobaric  centers  moving  westward  across  the  Antilles.  Riehl 
(1945)  provided  the  first  comprehensive  work  on  waves  in  the  easterlies,  stating  that 
these  waves  originated  over  Africa  as  a  result  of  barotropic  instability  in  the  easterly 
current.  He  found  that  these  waves  have  their  maximum  intensity  between  3,000  and 
4,500  m,  propagate  to  the  west  at  an  average  speed  of  5-8  ms-1,  and  have  a  period 
of  roughly  3-4  d.  The  top  of  the  moist  layer  to  the  west  of  the  trough  axis  can  be  as 
low  as  1,500  m  but  rises  rapidly  near  the  trough  axis,  attaining  a  maximum  height 
of  nearly  9,000  m  to  the  east  of  the  trough  axis  where  strong  convective  activity  is 
found. 

Riehl’s  original  work  was  expanded  upon  by  Merritt  (1964)  to  include  several 
variations  of  the  easterly  wave  model.  Also,  during  this  time,  Riehl’s  original 
model  began  to  draw  criticism.  Merritt  felt  this  controversy  resulted  from  a  general 
tendency  to  label  any  disturbance  in  the  easterly  flow  as  an  easterly  wave.  Atkinson 
(1971)  agreed  with  this  assessment,  claiming  that  many  forecasters  used  the  easterly 
wave  as  a  “crutch”  to  attempt  to  explain  all  tropical  weather  occurrences.  Further 
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controversy  arose  concerning  the  development  of  tropical  cyclones  from  waves  in  the 
easterlies.  While  Riehl  maintained  that  tropical  storm  development  resulted  from 
the  intensification  of  waves  in  the  easterlies,  Sadler  (1967)  led  the  opposing  view 
stating  that  all  tropical  vortices  form  initially  in  a  shear  zone  between  two  currents  of 
opposite  direction.  Despite  the  controversy  concerning  an  exact  definition  for  waves 
in  the  easterlies  and  their  role  in  tropical  cyclone  formation,  there  is  little  doubt  that 
perturbations  in  the  easterly  flow  do  exist  and  can  produce  significant  variations  in 
convective  activity  throughout  the  Caribbean  region. 

The  use  of  spectral  analysis  techniques  during  the  late  1960s  led  to  the  discovery 
of  atmospheric  oscillations  with  periods  longer  than  10  d.  Wallace  and  Chang  (1969) 
used  these  spectral  analysis  techniques  to  further  study  the  characteristics  of  waves  in 
the  easterlies.  Their  findings  confirmed  the  existence  of  fluctuations  in  the  zonal  and 
meridional  wind  components  with  a  periodicity  of  approximately  4-5  d  for  several 
tropical  locations,  including  Balboa,  Panama.  Similar  fluctuations  were  noted  for 
relative  humidity  and  surface  pressure.  An  additional  interesting  result  was  the 
detection  of  a  longer-period  oscillation  in  these  variables  having  a  period  greater 
than  10  d,  which  was  strongest  at  western  Pacific  stations.  This  oscillation  was 
downplayed  by  the  authors  and  was  attributed  to  the  sensitivity  of  the  spectra  in  the 
low-frequency  end  to  minor  changes  in  the  sampling  period.  It  was  felt  that  if  the 
period  of  study  had  been  changed  by  a  month  or  two  the  spectral  peak  might  have 
been  shifted  or  completely  eliminated. 

The  first  real  evidence  of  a  tropical  oscillation  with  a  period  significantly  longer 
than  a  wave  in  the  easterlies  was  presented  by  Madden  and  Julian  (1971).  Their 
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spectral  analysis  of  daily  rawinsonde  data  for  Canton  Island  (3°S,  172°W)  revealed 
a  41-53  d  fluctuation  in  the  85  and  15  kPa  zonal  winds,  tropospheric  temperature, 
and  station  pressure  in  the  tropical  central  Pacific  region.  Madden  and  Julian  (1972) 
expanded  their  original  study  to  include  several  tropical  stations  around  the  globe 
in  order  to  determine  the  spatial  extent  of  this  oscillation.  Their  results  revealed  an 
eastward-propagating,  large-scale  (zonal  wavenumber  one  or  two)  oscillation  with  a 
period  of  roughly  40-50  d.  This  oscillation  was  evident  in  station  pressure  and  zonal 
winds  for  all  stations  within  10°  latitude  of  the  equator  from  Singapore  eastward  to 
Curacao.  Balboa,  Panama  was  included  in  this  study  and  showed  spectral  peaks  in 
pressure  and  upper  tropospheric  zonal  winds  at  or  near  40-50  d. 

As  the  use  of  satellite  data  increased  during  the  1970s,  cloud  cover  and  outgoing 
longwave  radiation  (OLR)  data  began  to  be  used  to  detect  low-frequency  oscillations 
in  the  tropics.  Using  satellite  cloud  data  and  geopotential  height,  Yasunari  (1980) 
found  that  intraseasonal  variation  of  the  Indian  monsoon  during  northern  summer 
is  dominated  by  the  40-50  d  oscillation.  Murakami  (1980)  found  significant  15-30  d 
peaks  in  OLR  data  over  regions  of  southeast  Asia  during  winter,  while  Lyons  (1981) 
found  similar  spectral  peaks  over  Indonesia  and  the  equatorial  Pacific  regions.  These 
studies  demonstrated  the  ability  of  OLR  data  to  represent  convective  activity  in  the 
tropics,  with  small  OLR  values  indicating  areas  of  high  cloudiness  (associated  with 
convective  activity)  and  larger  OLR  values  showing  relatively  cloud-free  areas.  The 
success  of  those  studies  led  to  a  series  of  papers  which  used  OLR  data  to  detect  low- 
frequency  tropical  oscillations,  including  studies  by  Weickmann  (1983);  Weickmann, 
Lussky,  and  Kutzbach  (1985);  Lau  and  Chan  (1983);  and  Knutson,  Weickmann,  and 


6 


Kutzbach  (1986).  The  results  of  these  studies  confirmed  the  original  description  of 
the  oscillation  as  stated  by  Madden  and  Julian  and  showed  that  the  OLR  anomalies 
have  a  characteristic  propagation  speed  of  3-6  ms-1  and  are  strongest  over  the  tropical 
Indian  and  the  western  Pacific  oceans. 

While  most  of  the  studies  mentioned  above  have  concentrated  on  low-frequency 
oscillations  over  the  Indian  Ocean  and  western  Pacific  region,  Knutson  and  Weick- 
mann  (1987)  state  that  some  potentially-significant  OLR  fluctuations  seem  to  be 
occurring  over  parts  of  Central  and  South  America.  To  date,  very  little  research  has 
dealt  with  the  detection  and  significance  of  low-frequency  oscillations  in  this  tropical 
region.  They  also  go  on  to  state  that  the  practical  importance  of  the  oscillation, 
in  terms  of  relating  the  low-frequency  oscillations  to  a  local  forecast,  is  a  topic  that 
needs  to  be  investigated.  This  research  will  address  both  subjects. 

The  majority  of  work  done  in  analyzing  tropical  oscillations  to  date  has  focused 
on  atmospheric  variables  such  as  OLR,  geopotential  height,  temperature,  and  zonal 
and  meridional  wind  measurements.  Hartmann  and  Gross  (1988)  were  the  first 
to  detect  low-frequency  oscillations  using  daily  precipitation  records.  Their  results 
showed  significant  spectral  peaks  in  the  49-50  d  range  for  several  Indonesian  and 
western  Pacific  locations.  This  study  was  also  significant  in  that  it  began  to  focus 
attention  on  the  effect  these  atmospheric  oscillations  can  have  on  day-to-day  weather 
conditions  such  as  daily  rainfall  amounts  at  specific  locations.  Once  it's  been 
determined  that  these  oscillations  can  significantly  affect  rainfall  and  cloud  cover  at  a 
location,  the  next  important  step  is  to  develop  techniques  to  forecast  the  day-to-day 
changes  associated  with  these  oscillations. 
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CHAPTER  III 

RESEARCH  OBJECTIVES 

The  main  research  objective  was  to  examine  convective  variability  over  Panama 
through  spectral  analysis  of  OLR  time  series  data.  Any  significant  periodicities  in  the 
OLR  time  series  would  reflect  periodic  variations  of  cloud  cover  and  convection  over 
the  Panama  region.  The  spatial  extent  of  this  oscillation  would  then  be  determined 
through  the  construction  of  composite  OLR  anomaly  diagrams  over  a  larger  domain 
(120°W-40°W,  35°S-35°N)  corresponding  to  four  phases  (maximum,  minimum,  and 
two  inflection  points)  of  the  oscillation  cycle.  By  following  the  four  phases  of  the 
oscillation,  a  direction  of  movement  of  the  OLR  anomalies  could  be  inferred.  Similar 
composites  of  rawinsonde  data  corresponding  to  the  four  parts  of  the  oscillation  would 
be  constructed  which,  along  with  satellite  imagery,  would  allow  for  a  rough  analysis 
of  the  associated  synoptic  pattern. 

Another  objective  was  to  identify,  quantitatively  (through  use  of  rawinsonde 
wind  and  geopotential  height  composites)  and  qualitatively  (through  use  of  satellite 
imagery),  the  atmospheric  conditions  that  would  forecast  any  change  in  cloud  cover 
and  convective  activity  over  Panama.  By  maintaining  daily  continuity  on  these 
conditions,  operational  weather  forecasters  in  Panama  would  be  able  to  distinguish 
between  each  phase  of  the  oscillation,  allowing  them  to  significantly  improve  their 
medium-  to  long-range  forecasts. 

An  additional  goal  was  to  examine  the  correlation  between  daily- averaged  OLR 
values  over  Panama  and  daily  rainfall  totals  from  the  dense  network  of  raingauges  in 
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the  Panama  Canal  region.  A  comparison  of  the  rainfall  totals  and  OLR  values  would 
help  determine  the  usefulness  of  OLR  values  in  estimating  rainfall  over  a  region. 


CHAPTER  IV 


DATA 


OLR 

A  2-y  subset  (1  January  1984  through  31  December  1985)  of  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA)  Outgoing  Longwave  Radiation 
(OLR)  data  set  provided  the  primary  tool  used  in  this  research.  This  data  set  was 
established  in  June  1974  and  since  that  time  numerous  studies  have  used  these  OLR 
values  to  study  convective  activity  in  the  tropics.  These  OLR  data  were  derived  from 
window  radiance  measurements  obtained  from  a  scanning  radiometer  flown  aboard  the 
NOAA  7  polar-orbiting  satellite.  The  NOAA  7  is  one  in  a  series  of  sun- synchronous 
satellites  and  has  equator  crossing  times  of  0230  and  1430  LST.  The  spatial  resolution, 
as  originally  measured,  was  8  km  but  the  observations  have  been  spatially-averaged, 
digitized,  and  stored  in  a  2.5°  latitude-longitude  array.  A  more  detailed  description 
of  the  process  is  provided  by  Gruber  and  Winston  (1978). 

A  brief  summary  of  the  procedures  used  to  obtain  the  OLR  values  follows: 

The  Temperature  Humidity  Infrared  Radiometer  (THIR)  on  the  NOAA  7 
satellite  senses  window  radiance  (Ea)  in  the  infrared  window  region  of  11.5  to  12.5  /xm. 
Radiance  equivalent  brightness  temperatures  ( Trb )  are  obtained  by  solving  Planck’s 
law 

E  ClX^ 

*  C2/eT'i>*  —  1 


(1) 
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for  Trf, .  The  result  is 

rp  _  _ C2^ _ 

r6  ~  ln(ClA-5£^ 1  +  1)  ( 

where  A  is  the  wavelength  in  //m,  ci  =  3.74  x  10  16Wm2,  and  c 2  =  1.44  x  10  2mK. 

At  this  point,  an  empirical  regression  formula  developed  by  Ohring,  Gruber, 
and  Ellingson  (1984)  and  refined  by  Gruber  and  Krueger  (1984)  is  used  to  convert 
Trf,  to  flux  equivalent  brightness  temperatures  ( Tfb ): 

Tfb  =  aTx  +  bTl  (3) 

where  Tx  =  -3.66K  +  1.015(Trfc),  a  =  1.249,  and  b  =  0.001055K-1. 

Finally,  through  use  of  the  Stefan- Boltzman  law 


F0lr  —  <7  T ^ 


(4) 


where  a  =  5.67  x  10-8Wm-2K-4,  the  actual  OLR  flux  emittance  values  ( F0[r )  used 
in  this  study  were  obtained. 

A  relationship  between  OLR  values  and  equivalent  cloudtop  temperatures  ( Te ) 
provided  by  Cahalan,  Short,  and  North  (1982)  states 


AOLR 

ATe 


2Wm-2K-1 


with  a  reference  point  of  OLR  %  250  Wm'2  at  Te  =  273  K. 

The  assumption  made  in  using  these  data  is  that  small  values  of  OLR  represent 
areas  of  deep  convection  and  cloud  cover,  while  large  OLR  values  represent  relatively 
clear  conditions.  This  leads  to  the  major  disadvantage  of  using  OLR  data,  which 
is  that  it  is  only  one- dimensional  and  does  not  allow  for  an  analysis  of  the  vertical 
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structure  of  the  system.  For  example,  a  small  OLR  value  of  130  Wm~2  would  suggest 
that  convective  activity  and  precipitation  were  occurring  at  a  location.  However,  that 
low  OLR  value  could  be  reflecting  only  the  cold  cloud  tops  of  a  cirrus  anvil  associated 
with  convective  activity  occurring  far  from  that  location. 

Although  the  objective  of  this  study  was  to  examine  convective  variability  over 
Panama,  it  became  apparent  during  the  research  that  the  factors  affecting  convection 
in  that  area  were  originating  well  outside  of  the  Panama  region.  The  original  area 
of  study  continued  to  expand,  with  the  eventual  domain  of  the  OLR  field  shown  in 
Fig.  1. 

Upper- Air 

The  upper-air  data  used  in  this  study  are  from  Albrook  Air  Force  Station, 
Panama  (see  Fig.  2)  for  the  time  period  1  January  through  31  December  1984.  Vertical 
soundings  were  taken  twice  daily  at  the  site,  at  0000  and  1200  GMT.  Only  the  1200 
GMT  soundings  are  used  in  this  study  since  thunderstorms  are  often  still  in  the  local 
area  during  the  0000  GMT  run  and  can  affect  the  ability  of  the  sounding  to  represent 
conditions  on  a  larger  spatial  scale. 

Three  variables  were  used  from  the  soundings:  u-  and  v-wind  components  and 
height  of  the  significant  pressure  surfaces.  The  wind  data  were  obtained  from  pilot- 
balloon  (PIBAL)  observations  while  the  height  data  were  obtained  from  radiosonde 
observations  (RAOBS).  Missing  data  for  any  of  the  variables  were  resolved  through 


linear  interpolation  over  time. 


Figure  2.  Map  of  Panama  (adapted  from  Map  of  Panama,  U.  S.  Central  Intelligence  Agency) 
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Rainfall 

Rainfall  data  for  1984  and  1985  were  provided  by  the  Meteorology  and  Hy¬ 
drology  Branch  of  the  Panama  Canal  Commission.  Their  network  of  26  automated 
raingauges  provides  a  complete  coverage  of  rainfall  within  the  Panama  Canal  region. 
The  exact  locations  of  these  raingauges  are  shown  in  Fig.  3.  For  this  study,  daily 
rainfall  totals  from  each  raingauge  site  were  summed  to  provide  one  daily  rainfall 
total  for  the  entire  Canal  region.  Analysis  of  this  time  series  of  daily  rainfall  totals  is 
discussed  in  Chapter  VI. 
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CHAPTER  V 

DETERMINATION  OF  PERIODIC 
MODES  OF  CONVECTIVE  ACTIVITY 

Selection  of  Area  of  Study 

Periodic  fluctuations  in  an  area’s  OLR  time  series  are  easier  to  detect  if  an  area 
has  relatively  small  average  annual  OLR  values,  indicating  that  large  amounts  of  cloud 
cover  are  present  throughout  the  year.  Average  OLR  values  of  the  domain  for  the 
period  1  January  1984  through  31  December  1985  are  shown  in  Fig.  4.  These  average 
OLR  values  were  calculated  for  each  2.5°  latitude-longitude  grid  point  using  both 
the  0230  and  1430  LST  observations.  Note  the  relative  OLR  minima  extending  from 
central  Mexico  southeastward  through  Central  America.  The  lowest  OLR  values  in 
all  of  Central  America,  as  well  as  some  of  the  lowest  OLR  values  in  the  entire  domain, 
are  found  in  the  Panama  area,  especially  in  the  eastern  portions  of  the  country.  Also 
note  the  relativ<iijj  low  OLR  values  associated  with  the  trade  wind/monsoon  trough 
beginning  at  10°N,  120°W  and  extending  eastward  into  extreme  eastern  portions  of 
Colombia. 

In  addition  to  having  low  average  OLR  values,  an  important  criterion  in 
selecting  an  area  of  study  is  the  total  variance  of  OLR.  If  a  region  is  cloudy  all 
year  but  has  little  or  no  fluctuations  in  the  cloud  cover,  an  analysis  of  that  OLR  time 
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series  would  not  detect  any  significant  periodicities.  The  variance  was  calculated  for 
each  grid  point  in  the  domain  for  the  same  2-y  period  using  the  formula 

n 

»J  =  £(yi-F>!/( »-!)  (5) 

1=1 

where  Y{  represents  the  daily  averages  of  the  0230  and  1430  LST  observations,  Y  is  the 
2-y  average  OLR  for  that  point,  and  n  equals  the  total  number  of  daily  average  OLR 
values  over  the  2-y  period,  which  was  731.  In  these  calculations  the  average  of  the  two 
daily  OLR  observations  was  used  to  eliminate  any  variance  resulting  from  the  diurnal 
cycle.  The  total  variance  for  the  domain  is  shown  in  Fig.  5.  It  can  been  seen  that 
the  highest  total  variance  in  the  entire  domain  is  a  2000  (Wm~2)2  contour  located 
just  to  the  southwest  of  Panama  City.  In  addition,  the  entire  country  of  Panama, 
as  well  as  the  surrounding  area,  has  some  of  the  highest  OLR  variance  values  in  the 
domain.  This  combination  of  low  average  OLR  values  and  maximum  OLR  variance 
should  make  the  Panama  area  daily  OLR  values  favorable  for  time  series  analysis. 

The  exact  boundaries  of  the  area  of  study  must  now  be  specified.  They  must  be 
broad  enough  to  be  representative  of  the  region  as  a  whole,  but  small  enough  so  that 
the  changes  in  the  OLR  field  that  could  be  associated  with  synoptic-scale  features 
are  not  averaged  out.  Figure  6  shows  the  portion  of  the  domain  selected  for  analysis, 
with  the  boundaries  extending  from  5°-12.5°N  and  85°-75°W.  These  boundaries  were 
selected  on  the  basis  of  the  average  OLR  values  and  total  OLR  variance  as  shown  in 
Figs.  4  and  5,  as  well  as  the  fact  that  this  area  is  subject  to  the  influence  of  large- 
scale  effects  such  as  northward  migrations  of  the  monsoon  trough,  the  passage  of 
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waves  in  the  easterlies,  and  the  occasional  southward  advance  of  shear  lines  towards 
the  Panama  area. 

Methods  of  Analysis 

There  are  20  OLR  grid  points  contained  in  the  area  of  study  shown  in  Fig.  6. 
A  single  daily  OLR  value  was  calculated  for  the  area  of  study  by  taking  the  average 
of  the  20  0230  LST  and  the  20  1430  LST  OLR  point  values.  A  study  of  the  area- 
averaged  OLR  as  opposed  to  a  single  grid  point  greatly  reduces  the  noise  in  the  OLR 
time  series,  allowing  a  differentiation  between  local  and  large-scale  effects.  These 
area-averaged  daily  OLR  values  comprise  the  OLR  time  series  for  the  Panama  area 
and  will  be  used  in  the  time  series  analysis.  Plots  of  these  time  series  of  daily  OLR 
averages  for  1984  and  1985  are  shown  in  Figs.  7  and  8. 

Division  of  Data  into  Dry  and  Wet  Seasons 

The  next  step  in  preparing  the  data  for  analysis  is  to  recognize  that  the 
conditions  in  the  Panama  area  are  not  constant  throughout  the  year  but  are  divided 
into  dry  and  wet  seasons.  Several  studies  over  the  past  20  y  have  concerned  themselves 
with  detecting  periodicities  in  tropical  variables  and  have  chosen  not  to  deal  with 
these  two  seasons  separately.  This  would  not  be  an  appropriate  method  of  study 
for  the  Panama  area  since  the  two  seasons  are  so  distinctly  different  in  terms  of 
cloud  cover  and  precipitation  as  well  as  the  upper-level  wind  patterns  associated  with 
these  features.  A  periodicity  that  is  detected  in  both  the  dry  and  wet  seasons  could 
be  the  reflection  of  two  completely  different  synoptic  or  local  conditions  depending 
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on  the  time  of  the  year.  For  this  reason,  the  2-y  period  has  been  divided  into  dry 
and  wet  seasons  through  the  subjective  analysis  of  raw  OLR  values  as  well  as  daily 
precipitation  totals.  Both  time  series  show  distinct  delineations  between  the  two 
seasons,  which  are  as  follows: 

•  Dry  season  1984:  1  January-9  May  (130  days) 

•  Wet  season  1984:  10  May-4  December  (209  days) 

•  Dry  season  1985:  5  December  1984-2  May  1985  (149  days) 

•  Wet  season  1985:  3  May-8  December  (220  days). 

Removal  of  Seasonal  Trends 

Although  a  distinct  discontinuity  in  the  raw  OLR  and  precipitation  data  marks 
each  change  of  season,  a  long-period  seasonal  trend  is  clearly  visible  in  each  season. 
This  can  be  seen  in  the  OLR  time  series  shown  in  Figs.  7  and  8.  In  preparing 
each  season’s  OLR  time  series  for  spectral  analysis,  Chatfield  (1980)  recommends 
removing  any  obvious  long-period  seasonal  trends  from  the  time  series  so  that  this 
trend  does  not  detract  from  the  power  of  any  higher-frequency  oscillations.  To 
eliminate  these  seasonal  trends,  a  2nd  degree  polynomial  trend  remover  method  was 
used,  as  described  by  Neiswanger  (1948).  This  method  fits  the  data  to  a  polynomial 
curve,  with  the  fitted  function  providing  a  measure  of  the  trend.  The  residuals  about 
this  fitted  curve  provide  an  estimate  of  local  fluctuations  which  are  independent  from 
the  long-period  trend.  These  fluctuations  about  the  fitted  curve  (which  has  been 
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set  to  zero)  for  each  season’s  OLR  time  series  are  shown  in  Figs.  9,  10,  11,  and  12. 
Note  that  the  broad  curves  associated  with  seasonal  trends  in  the  OLR  time  series 
of  Figs.  7  and  8  have  been  removed  and  the  daily  OLR  values  are  now  represented 
as  deviations  about  the  fitted  curves.  These  four  OLR  time  series  are  now  ready  for 
time  series  analysis. 


Time  Series  Analysis 

The  method  used  to  determine  whether  any  significant  periodicities  were  present 
in  the  OLR  time  series  was  harmonic  analysis,  which  is  the  process  of  decomposing 
a  time  record  into  the  sum  of  sinusoids  of  various  frequencies  and  amplitudes.  The 
idea  behind  harmonic  analysis  is  that  a  complicated  time  series  can  be  understood  by 
break.  ig  it  up  into  these  simple  frequency  components  and  studying  them  separately. 
A  detailed  discussion  of  harmonic  analysis  is  provided  by  Panofsky  and  Brier  (1968). 
One  type  of  harmonic  analysis  which  is  used  to  decompose  this  OLR  time  series  into 
its  sine  and  cosine  components  is  the  fast  Fourier  transform  (FFT)  method.  The 
International  Mathematical  and  Statistical  Library  (IMSL)  User’s  Manual  describes 
the  FFT  method  and  provides  a  computer  routine  that  uses  the  FFT  method  to 
perform  a  periodogram  analysis.  This  periodogram  analysis  calculates  the  sum  of 
the  squared  amplitudes  of  the  sinusoids,  or  spectral  density,  for  each  frequency.  The 
IMSL  periodogram  analysis  was  done  for  the  box-averaged  OLR  time  series  of  each 
season,  with  the  results  shown  in  Figs.  13  and  14. 

As  seen  in  Fig.  13,  dry  season  1984  shows  a  spectral  peak  at  around  65  d, 
but  the  largest  spectral  peak  for  this  season  is  concentrated  in  the  11-12  d  range. 


i - 1  i - 1 - i  i - 1 - 1  i - 1  i  :  r 

09  09  0*  0E  02  01  0  01-  02-  OS-  Ofr-  09-  09- 

3*.IN/M  U10 


30 


PERXOOOGRAM 
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Figure  13 


Numerical  output  from  IMSL  periodogram  analysis  for  dry  season  and 
wet  season  1984.  Periods  are  in  days.  I(w(k))  represents  spectral  density. 
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PERIODOGRAM  ANALYSIS  -  DRY  SEASON  1985  PERIODOGRAM  ANALYSIS  -  WET  SEASON  1985 
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Numerical  output  from  IMSL  periodogram  analysis  for  dry  season  and 
wet  season  1985.  Periods  are  in  days,  I(w(k))  represents  spectral  density. 
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On  all  of  the  analyses  the  spectral  density  begins  to  decrease  rapidly  below  6  d,  so 
those  portions  of  the  analyses  have  been  omitted.  The  periodogram  analysis  for  wet 
season  1984  shows  a  spectral  peak  of  approximately  35-40  d,  which  was  not  evident 
in  the  dry  season.  A  spectral  peak  is  also  found  at  6-7  d,  however,  another  high 
concentration  of  spectral  density  in  this  season  can  again  be  found  at  around  12  d. 
The  largest  spectral  peak  in  dry  season  1985  is  found  at  approximately  25  d,  but 
another  large  peak  is  evident  at  11-12  d.  Finally,  wet  season  1985  shows  a  spectral 
peak  in  the  60-70  d  range,  but  the  strongest  peak  again  appears  to  be  located  at 
roughly  11-12  d. 

The  periodogram  analyses  of  Figs.  13  and  14  show  there  are  several  frequencies 
that  display  spectral  peaks  through  the  2-y  period,  but  by  far  the  strongest  and  most 
consistent  is  found  in  the  12-d  range.  The  next  step  is  to  apply  a  different  analysis 
technique  to  the  data  to  see  if  the  results  are  consistent.  Several  additional  techniques, 
as  well  as  a  comprehensive  discussion  on  time  series  analysis,  are  presented  by  Newton 
(1988).  One  of  the  methods  outlined  in  his  book  TIMESLAB  is  the  correlogram.  The 
correlogram  is  a  plot  of  the  autocorrelation  coefficient  (p(i/))  vs  the  lag  period  (i/). 
In  other  words,  it  is  measuring  the  correlation  of  the  data  with  themselves  except 
“lagged”  by  a  certain  number  of  time  units.  This  technique  was  applied  to  the  OLR 
time  series  for  each  season  and  the  results  are  shown  in  Figs.  15  and  16.  Although  the 
p( v)  trace  never  exceeds  the  95%  confidence  band,  the  correlogram  for  dry  seasons 
1984  and  1985  exhibit  sinusoidal  patterns  with  a  periodicity  in  the  time  frame  of 
around  12  d.  The  correlograms  for  both  wet  seasons  indicate  a  weak  12-d  periodicity, 
but  not  nearly  as  strong  as  the  dry  seasons. 
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DRY  SEASON  1984 


68  12.66  24.66  36.66  48.66  66.66 


WET  SEASON  1984 


I .  00' - i - 1 - i - i - 1 - 1 - 1 - t - 1 - ‘ 

6.66  12.66  24.66  36.66  49.60  66.60 


Figure  15.  OLR  correlograms  and  95%  confidence  bands  for  1984.  Y  axes  represent 
autocorrelation  coefficient  ( p(i /)),  X  axes  represent  lag  period  in  days  (u). 
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DRY  SEASON  1985 


0.00  12.00  24.00  36.00  48.00  60.00 


WET  SEASON  1985 


0.00  12.00  24.00  36.00  48.00  60.00 


Figure  16.  OLR  correlograms  and  95%  confidence  bands  for  1985.  Y  axes  represent 
autocorrelation  coefficient  (p(«/)),  X  axes  represent  lag  period  in  days  ( i /). 
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In  summary,  the  two  time  series  analysis  techniques  used  to  determine  any 
periodicities  in  the  OLR  time  series  were  the  IMSL  periodogram  and  the  TIMESLAB 
correlogram.  Strong  spectral  power  in  the  12-d  period  was  evident  in  all  four  seasons 
in  the  periodogram  analysis,  while  the  correlogram  showed  strong  12-d  fluctuations  in 
the  dry  seasons  of  both  years  with  weaker  12-d  fluctuations  evident  in  the  wet  seasons. 
Other  spectral  peaks  were  present  during  some  seasons  but  were  not  as  strong  or  as 
persistent  as  the  12-d  fluctuation. 
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CHAPTER  VI 

ANALYSIS  OF  DOMINANT  MODE  OF 
CONVECTIVE  ACTIVITY  OVER  PANAMA 

Examination  of  12-Day  Oscillation 

Using  the  spectral  analysis  techniques  described  in  Chapter  V,  the  daily- 
averaged  OLR  time  series  over  the  Panama  region  have  all  shown  spectral  peaks 
of  approximately  12  d.  To  determine  the  actual  strength  of  this  12-d  oscillation,  a 
band-pass  recursive  filter  developed  by  Shanks  (1967)  and  used  by  Murakami  (1980) 
was  applied  to  the  OLR  time  series  of  Figs.  9,  10,  11,  and  12.  This  recursive  filter 
computes  each  output  point  as  a  weighted  sum  of  input  points  plus  a  weighted  sum 
of  previously  computed  output  points.  This  filter  allows  us  to  set  the  periods  where 
full  (1.0)  and  lower  (0.5)  amplitude  responses  are  desired.  In  this  study  the  two 
periods  chosen  were  12  d  (full  response)  and  9  d  (0.5  response),  which  translate  into 
frequencies  (measured  by  27r/period)  wo  and  w i  respectively. 

The  filtered  time  series  data  at  day  n  (fn)  were  obtained  from  the  unfiltered 
data,  fn,  by 

fn  =  a(fn  fn  — 2)  ^lfn—  1  ^fn— 2» 

n  =  l,2,...,iV,  (6) 

4  +  2A  n  +  nl' 


where 
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bi  = 


b2  = 


2(0} -4) 

4  +  2AO  + 

4  -  2Aft  +  n20 
4  +  2AQ  +  &' 


(7) 


and 


Afi  =  2 


sinwi 


sinu>2 


nl  = 


1  +  coswi  1  +  cosu>2 

4(sinu>i)(sinu>2) 

(1  +  coswi)(l  +  cosu>2) 


(8) 


The  frequencies  aq  and  uj2  in  Eq.  8  are  the  frequencies  at  which  we  get  0.5 
response  on  both  sides  of  wo-  Recall  we  have  already  specified  the  1.0  response 
frequency  for  12  d  (wo),  and  the  0.5  response  frequency  for  9  d  (wi).  To  obtain  the 
other  0.5  response  frequency  (u)2)  the  following  relationship  is  used: 


u>2  =  Wq/wi. 


(9) 


The  filtered  data  are  obtained  by  first  applying  Eq.  (6)  to  the  original  time 
series  data  (from  which  the  mean  and  any  linear  trends  have  been  removed)  from 
n  =  3  to  N  by  assuming  that  fj  =  f2  =  0  to  obtain  tentative  mean  values  for  fn. 
Then,  using  these  tentative  mean  values,  the  process  is  reversed  in  time  (from  N  —  2 
to  n  =  3)  to  obtain  the  final  fn  values.  A  significant  advantage  to  using  this  procedure 
is  that  it  results  in  a  zero  phase  shift  for  all  frequencies. 

The  amplitude  response  (Rn)  for  each  frequency  is  given  by 


Rn  =  a(l  -  *})/(!  +  b\zn  +  b2zl), 


(10) 
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where  zn  =  e~^2r^n\n  =  the  period  in  days,  and  a, 61,  and  62  are  calculated  from 
Eq.  (7).  The  actual  amplitude  response  curve  of  the  12-d  filter  used  in  this  study  for 
a  length  of  n  =  100  d  is  shown  in  Fig.  17. 

The  plots  of  the  12-d  filtered  OLR  time  series  for  each  season  are  shown  in 
Figs.  18,  19,  20,  and  21.  Note  that  in  all  seasons  significant  amounts  of  the  OLR 
deviations  have  been  retained  after  filtering  for  12  d.  A  quantitative  measure  of  the 
strength  of  the  signal  can  be  found  by  comparing  the  standard  deviations  of  the 
unfiltered  and  filtered  time  series.  The  variances  for  the  unfiltered  OLR  time  series 
of  Figs.  9,  10,  11,  and  12  have  been  calculated  as  334.9,  497.3,  156.3,  and  449.4 
(Wm-2)2,  respectively.  In  comparison,  the  12-d  filtered  OLR  time  series  of  Figs.  18, 
19,  20,  and  21  have  variances  of  98.0,  118.8,  29.2,  and  94.1  (Wm-2)2,  so  roughly  20- 
25%  of  the  variance  of  the  original  OLR  time  series  has  been  retained  after  filtering 
for  12  d.  An  examination  of  these  diagrams  shows  that  while  this  12-d  oscillation  is 
evident  in  each  season,  it  is  strongest  and  most  consistent  during  the  1984  dry  and 
wet  seasons.  Of  particular  interest  is  the  1984  wet  season  filtered  time  series,  which 
shows  an  apparent  modulation  of  the  12-d  oscillation  occurring  every  70-80  d.  This 
feature  will  be  addressed  in  Chapter  VII. 

OLR  Composite  Diagrams 

The  majority  of  this  research  to  this  point  has  concentrated  on  detecting 
periodicities  in  the  OLR  field  over  the  Panama  region.  Now  that  a  12-d  oscillation 
has  been  found,  the  next  step  is  to  determine  the  spatial  extent  of  that  oscillation. 
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This  was  accomplished  through  the  use  of  OLR  composite  diagrams,  which  were 
constructed  as  follows: 

•  The  OLR  time  series  of  each  2.5°  OLR  grid  point  in  the  domain  of  Fig.  1  (957 
total  points)  was  divided  into  the  same  dry  and  wet  seasons  as  the  Panama  area 
for  1984  and  1985. 

•  The  seasonal  mean  and  long-period  trends  were  removed  from  the  OLR  seasonal 
time  series  at  each  point. 

•  The  OLR  time  series  for  each  point  was  filtered  for  12  d  in  the  same  manner  as 
the  Panama  area  OLR  time  series. 

•  Cutoffs  of  ±1  standard  deviation  were  established  for  each  of  the  12-d  filtered 
OLR  time  series  shown  in  Figs.  18,  19,  20,  and  21.  For  each  season,  any  day  in 
which  the  OLR  deviation  in  the  Panama  area  time  series  exceeded  one  standard 
deviation  above  the  seasonal  mean  was  considered  a  Maximum  OLR  day,  or  in 
other  words,  a  relatively  doud-free  day  over  the  Panama  area.  For  these  same 
days,  the  corresponding  OLR  deviations  from  the  12-d  filtered  OLR  time  series 
for  each  point  in  the  domain  were  collected  and  averaged.  This  produced  average 
OLR  deviations  for  each  point  in  the  domain  corresponding  to  Maximum  OLR 
days  in  Panama.  By  constructing  this  OLR  composite  diagram  we  can  now 
determine  whether  conditions  occurring  over  Panama  with  a  periodicity  of  12  d 
have  any  relationship  with  the  rest  of  the  domain.  The  same  procedure  was 
used  to  construct  composite  diagrams  corresponding  to  Minimum  OLR  days 
(OLR  deviations  less  than  one  standard  deviation  below  the  seasonal  mean) 


over  the  Panama  area. 
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•  Composite  diagrams  were  also  constructed  for  the  two  inflection  points  found 
in  the  12-d  filtered  time  series.  They  correspond  to  the  inflection  point  going 
from  cloud-free  to  cloudy  conditions  (Max-to-Min)  and  going  from  cloudy  to 
cloud-free  conditions  (Min-to-Max). 

As  was  mentioned  above,  each  individual  point  in  the  domain  was  filtered  for 
12  d  without  first  performing  a  spectral  analysis  to  determine  if  any  spectral  peaks 
existed  in  that  period.  The  reasoning  behind  this  procedure  is  that  if  no  12-d  signal 
existed  anywhere  else  in  the  domain,  application  of  a  12-d  filter  to  the  OLR  time 
series  at  each  point  would  filter  out  all  of  the  signal,  resulting  in  time  series  that 
would  show  little  or  no  OLR  deviations.  Any  composite  analysis  using  these  time 
series  would  result  in  a  diagram  that  would  also  show  very  small,  if  not  zero,  OLR 
anomalies.  As  will  be  seen  in  the  next  section,  this  was  not  the  case. 

Composite  diagrams  were  constructed  for  both  years,  but  since  the  12-d  oscil¬ 
lation  was  most  evident  in  the  1984  OLR  time  series,  only  that  year  will  be  shown 
here. 


Composite  Diagrams  for  Dry  Season  1984 
The  composite  diagrams  for  dry  season  1984  are  shown  in  Figs.  22,  23,  24, 
and  25.  Beginning  the  12-d  cycle  with  the  Maximum  OLR  composite  (Fig.  22),  a 
large  area  of  maximum  OLR  deviations  implying  relatively  clear  conditions  can  be 
seen  through  most  of  Central  America  extending  to  the  west,  along  a  latitude  band 
of  about  2°-3°N.  A  broad  area  of  negative  OLR  anomalies,  representing  more  cloud 
cover  than  usual,  is  centered  at  7°N,  105°W,  extending  beyond  the  domain  to  the 


Jim  composite  diagram  corresponding  to  Max -to- Mm  inflection  points  over  Panama  lor  dry  season 
ontoured  values  represent  deviations  from  seasonal  averages  in  units  of  Win  ‘ . 
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west.  Three  days  later,  the  Max-to-Min  inflection  point  composite  (Fig.  23)  shows 
the  positive  OLR  anomalies  over  the  Panama  area  have  diminished  and  appear  to 
have  shifted  slightly  to  the  north.  The  major  difference  in  this  composite  is  the 
development  of  the  strong  dipole  OLR  anomaly  centered  at  7°N,  115°W  and  5°S, 
105°W.  Note  the  intensity  of  the  minimum  values  in  the  northern  OLR  minimum 
region,  which  are  roughly  —  27Wm~2.  Moving  ahead  three  more  days  in  the  cycle  to 
the  Minimum  OLR  composite  (Fig.  24)  it  appears  that  these  strong  minimum  OLR 
anomalies  have  built  towards  the  east,  with  most  of  Central  America  now  dominated 
by  these  low  OLR  values.  The  remnants  of  the  strong  negative  OLR  anomalies  over 
the  eastern  Pacific  can  be  seen  extending  from  the  Central  American  region  to  the 
west  along  a  latitude  of  roughly  3°N.  At  around  7°N,  105°W  a  new  area  of  positive 
OLR  anomalies  can  be  seen  developing,  extending  to  the  west  beyond  the  boundary 
of  the  domain.  An  interesting  comparison  can  be  made  at  this  point  between  the 
Maximum  and  Minimum  composites.  It  appears  the  two  are  mirror  images  of  each 
other,  in  both  location  and  intensity  of  OLR  anomaly  centers.  This  comparison 
holds  for  the  inflection  point  composites  as  well,  as  strong  similarities  can  be  seen 
between  the  Min-to-Max  composite  and  the  Max-to-Min  composite.  The  Min-to-Max 
composite  (Fig.  25)  shows  the  large  negative  OLR  anomalies  over  Central  America 
have  dissipated  and  shifted  to  the  north,  while  a  large  dipole  positive  OLR  anomaly 
has  developed  in  the  western  portion  of  the  domain.  Again,  this  scenario  is  nearly 
the  exact  opposite  of  the  Max-to-Min  inflection  point  composite. 

In  summary,  it  appears  that  in  the  1984  dry  season  a  large-scale,  west-to- 
east  moving  feature  is  responsible  for  creating  large  OLR  anomalies  across  the 
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domain  which  extend  from  beyond  the  western  edge  of  the  domain  eastward  through 
Panama  and  Central  America.  Throughout  the  domain  these  anomalies  appear  to  be 
significant  in  both  size  and  strength,  suggesting  that  this  12-d  oscillation  is  affecting 
not  only  Panama,  but  much  of  the  eastern  Pacific  as  well. 

Composite  Diagrams  for  Wet  Season  1984 
Strong  OLR  anomalies  during  the  wet  season  of  1984  can  be  seen  on  the 
composite  diagrams,  however,  the  development  and  apparent  movement  of  these 
anomalies  appears  to  be  quite  different  from  the  dry  season.  Figures  26,  27,  28,  and 
29  correspond  to  OLR  composites  through  one  complete  cycle  of  the  12-d  oscillation. 
Beginning  with  the  Maximum  OLR  composite  (Fig.  26),  the  large  positive  OLR 
anomalies  can  be  seen  throughout  Central  America  extending  westward  to  the  edge  of 
the  domain.  An  area  of  low  OLR  values  is  centered  just  off  the  southwest  Mexico  coast 
at  15°N,  100°W,  extending  off  to  the  northeast.  This  differs  significantly  from  the 
negative  anomaly  pattern  associated  with  the  dry  season  maximum  OLR  composite 
diagram,  as  this  minimum  OLR  anomaly  is  located  farther  to  the  northeast.  More 
differences  in  the  OLR  anomaly  field  can  be  seen  three  days  later  in  the  Max-to-Min 
composite  diagram  (Fig.  27).  It  appears  the  large  positive  OLR  anomaly  over  Panama 
and  Central  America  has  kept  much  of  its  original  intensity  despite  a  slight  decrease  in 
size  and  has  shifted  to  the  west-northwest,  just  south  of  the  coasts  of  El  Salvador  and 
Guatemala.  The  negative  OLR  anomaly  that  was  located  off  the  southwest  Mexico 
coast  has  expanded  significantly  and  built  farther  to  the  west.  Again,  this  pattern 
is  quite  different  from  the  dry  season  pattern  in  which  the  OLR  anomalies  appeared 


omposite  diagram  corresponding  to  Maximum  OLR.  days  over  Panama  for  wet  season  1984.  Contoured 
represent  deviations  from  seasonal  averages  in  units  of  Wm 


Figure  29. 
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to  be  drifting  to  the  east.  Three  days  later,  the  Minimum  OLR  composite  (Fig.  28) 
shows  that  convection  has  developed  explosively  throughout  Central  American  with 
the  lowest  OLR  values  centered  over  the  Costa  Rica-Nicaragua  region.  Convection 
also  seems  to  have  intensified  along  the  entire  monsoon  trough.  Meanwhile,  the 
positive  OLR  anomaly  that  was  located  off  the  southern  coasts  of  El  Salvador  and 
Guatemala  has  diminished  somewhat  in  size  but  has  maintained  much  of  its  original 
intensity  and  has  shifted  to  the  west-northwest.  Finally,  the  Min-to-Max  inflection 
point  composite  (Fig.  29)  shows  that  the  large  area  of  negative  OLR  anomalies 
centered  over  Costa  Rica  and  Nicaragua  3  d  earlier  has  diminished  in  size  and  moved 
to  the  west-northwest.  One  feature  the  wet  season  pattern  has  in  common  with  the 
dry  season  is  that  the  Maximum  and  Minimum  composites,  as  well  as  the  Max-to- 
Min  and  Min-to-Max  inflection  point  composites,  appear  to  be  mirror  images  of  each 
other. 

In  comparison  with  the  dry  season  OLR  anomalies,  which  appeared  to  be 
originating  off  the  western  boundary  of  the  domain  and  building  from  west  to  east, 
the  source  region  for  the  wet  season  OLR  anomalies  appears  to  be  directly  over 
the  Panama  region  with  the  OLR  anomalies  drifting  slowly  to  the  west-northwest. 
However,  it  is  impossible  to  infer  a  possible  synoptic  scenario  associated  with  the 
development  and  movement  of  these  OLR  anomalies  without  first  examining  the 
vertical  structure  of  the  atmosphere.  Also,  by  examining  the  upper-air  sounding  for 
any  significant  changes  as  the  12-d  oscillation  moves  over  Panama,  forecasters  may  be 
able  to  use  these  upper-air  analyses  to  forecast  changes  in  cloud  cover  and  convection 


associated  with  the  oscillation. 
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Vertical  Structure  of  12-Day  Oscillation 

Upper-air  data  from  Albrook  Air  Force  Station,  Panama  were  studied  in  order 
to  determine  the  vertical  structure  of  the  12-d  oscillation,  as  well  as  to  determine  if 
the  movement  of  the  oscillation  through  Panama  could  be  detected.  There  were  two 
methods  used  to  obtain  these  goals:  First,  the  IMSL  periodogram  analysis  routine 
was  performed  on  all  three  variables  (u-  and  v-wind  component  and  height  of  pressure 
surfaces)  to  determine  if  any  spectral  power  in  the  12-d  range  was  present  for  any  of 
the  variables.  Second,  vertical  composites  for  each  of  the  four  phases  of  the  oscillation 
were  constructed  for  each  variable  to  determine  if  any  significant  deviations  from 
the  mean  were  evident.  If  these  deviations  were  consistent  over  the  period  of  the 
oscillation  they  could  be  used  as  a  forecasting  tool.  These  studies  were  divided  into 
dry  and  wet  seasons,  and  the  procedures  and  results  are  discussed  in  the  following 
sections. 

Dry  Season  1984  Vertical  Composites 

The  procedures  used  to  study  the  upper-air  data  were  similar  to  those  used  on 
the  OLR  time  series.  First,  u-  and  v-wind  components  and  height  of  the  pressure 
surfaces  were  obtained  from  10  levels  (100,  85,  70,  50,  40,  30,  25,  20,  15,  and  10  kPa). 
For  each  variable  at  every  level,  the  seasonal  mean  and  any  long-period  seasonal  trends 
were  removed,  and  the  IMSL  periodogram  analysis  routine  was  then  performed.  The 
results  from  the  analyses  (not  shown)  are  briefly  summarized  as  follows: 

•  Spectral  peaks  were  evident  at  all  levels  for  the  u-wind  component  in  the 
30-60  d  range,  consistent  with  the  findings  of  Madden  and  Julian  (1972). 
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Broad  spectral  peaks  at  15-30  d  were  also  evident  at  some  levels.  Some 
spectral  power  at  all  levels  was  present  in  the  11-12  d  range. 

•  Very  broad  spectral  peaks  in  the  30-60  d  range  were  present  at  some  levels 
for  the  v-wind  component,  with  another  strong  spectral  peak  detected  at 
around  20-25  d.  However,  the  strongest  and  most  persistent  period  was 
found  at  11-12  d. 

•  The  strongest  spectral  peaks  in  the  height  field  were  found  30-45  and  20- 
25  d  periods  in  the  upper  levels  (40-10  kPa).  Weak  spectral  power  was 
detected  at  all  levels  for  the  period  of  11-12  d. 

The  results  of  the  IMSL  analyses  have  shown  that  other  low-frequency  oscilla¬ 
tions  can  be  detected  in  each  of  the  three  variables  from  the  Panama  sounding.  The 
11-12  d  oscillation,  while  not  always  the  strongest  spectral  peak,  was  evident  at  each 
level  with  the  strongest  signal  found  in  the  v-wind  component. 

The  next  step  in  the  analysis  was  to  filter  each  variable  at  every  level  for  12  d. 
similar  to  the  process  used  on  the  OLR  data.  Once  this  was  accomplished,  these 
analyses  were  used  to  construct  vertical  composite  diagrams  corresponding  to  each  of 
the  four  phases  (Maximum,  Max-to-Min  inflection  point,  Minimum,  and  Min-to-Max 
inflection  point)  of  the  Panama  area  12-d  filtered  OLR  time  series  (Fig.  18).  These 
vertical  composites  would  show  the  deviation  from  the  seasonal  means  associated 
with  each  phase  of  the  12-d  oscillation.  Using  these  composites,  a  forecaster  should 
be  able  to  determine  which  phase  of  the  12-d  oscillation  Panama  was  in  simply  by 
maintaining  continuity  of  the  daily  1200  GMT  upper-air  analyses.  The  following  is  a 
brief  description  of  the  vertical  composites: 
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•  The  vertical  composites  for  the  u-wind  component  corresponding  to  each 
of  the  four  phases  of  the  oscillation  are  shown  in  Figs.  30,  31,  32,  and 
33.  The  largest  deviations  from  the  seasonal  average  on  the  composites 
are  only  around  2-3  kt,  but  more  important  is  the  trend  of  the  deviations 
through  the  cycle  for  the  levels  50  kPa  and  above.  Note  that  at  the 
Maximum  composite  (Fig.  30). the  winds  have  a  slightly  weaker  easterly 
component  from  50-30  kPa,  and  a  slightly  stronger  westerly  component 
above  30  kPa.  Moving  through  the  cycle  to  the  Minimum  composite 
(Fig.  32)  the  anomalies  have  reversed,  with  slightly  stronger  easterly  winds 
from  50-30  kPa  and  a  slightly  weaker  westerly  component  above  30  kPa. 
Although  this  only  amounts  to  a  reversal  in  the  wind  component  of  about 
4  kt,  the  fact  that  this  trend  can  be  detected  in  the  entire  upper-level 
u-wind  field  could  be  a  significant  feature  used  to  forecast  the  phase  of 
the  12-d  wave. 

•  The  v-wind  composites  are  shown  in  Figs.  34,  35,  36,  and  37.  A  stronger 
trend  can  be  seen  in  the  v-component  at  the  25-20  kPa  level.  In  the 
Maximum  composite  (Fig.  34)  the  v-wind  at  that  level  is  only  about  1  kt 
from  the  south,  which  is  roughly  5  kt  less  than  the  seasonal  average.  By 
the  time  the  12-d  oscillation  cycle  reaches  the  Minimum  phase  6  d  later 
(Fig.  36),  the  southerly  wind  component  has  increased  by  over  8  kt  at  the 
25-20  kPa  level. 

•  The  composites  of  the  height  of  the  pressure  surfaces  are  shown  in  Fig.  38. 
The  height  composites  differ  from  the  u-and  v-wind  composites  in  that 
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(EASTERLY)  KNOTS  (WESTERLY) 


Figure  30.  Vertical  composite  for  u-wind  component  corresponding  to  Maximum 
OLR  days,  dry  season  1984.  Solid  line  represents  seasonal  average, 
dashed  line  represents  deviation. 
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Figure  31.  Vertical  composite  for  u-wind  component  corresponding  to  Max- to- Min 
OLR  inflection  points,  dry  season  1984.  Solid  line  represents  seasonal 
average,  dashed  line  represents  deviation. 
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Figure  32.  Vertical  composite  for  u-wind  component  corresponding  to  Minimum 
OLR  days,  dry  season  1984.  Solid  line  represents  seasonal  average, 
dashed  line  represents  deviation. 
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(EASTERLY)  KNOTS  (WESTERLY) 


Figure  33.  Vertical  composite  for  u-wind  componeni  corresponding  to  Min-to-Max 
OLR  inflection  points,  dry  season  1984.  Solid  line  represents  seasonal 
average,  dashed  line  represents  deviation. 
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Figure  34.  Vertical  composite  for  v-wind  component  corresponding  to  Maximum 
OLR  days,  dry  season  1984.  Solid  line  represents  seasonal  average, 
dashed  line  represents  deviation. 
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(NORTHERLY)  KNOTS  (SOUTHERLY) 


Figure  36.  Vertical  composite  for  v-wind  component  corresponding  to  Minimum 
OLR  days,  dry  season  1984.  Solid  line  represents  seasonal  average, 
dashed  line  represents  deviation. 
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Figure  37.  Vertical  composite  for  v-wind  component  corresponding  to  Min-to-Max 
OLR  inflection  points,  dry  season  1984.  Solid  line  represents  seasonal 
average,  dashed  line  represents  deviation. 


HEIGHT  DEV  (METERS) 


Figure  38.  Vertical  composite  of  height  anomalies  for  dry  season  1984.  Solid  line 
along  zero  deviation  line  represents  seasonal  mean,  solid  traces  represent 
deviations  from  the  mean.  Trace  1  corresponds  to  Maximum  OLR 
composite,  Trace  2  to  Max-to-Min  OLR  composite,  Trace  3  to  Minimum 
OLR  composite.  Trace  4  to  Min-to-Max  OLR  composite. 
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the  height  deviations  can  be  detected  at  all  levels,  as  opposed  to  the 
wind  composites  which  were  evident  at  only  the  upper  levels.  Beginning 
with  the  Maximum  composite  (Trace  1),  the  upper  levels  show  maximum 
height  anomalies  while  the  lower  levels  remain  near  the  seasonal  mean. 
Three  days  later  (Trace  2),  the  upper  levels  have  shifted  towards  the 
seasonal  mean  while  the  lower  levels  show  a  relatively  large  negative  height 
anomaly.  The  minimum  composite  (Trace  3)  shows  that  while  the  lower 
levels  retreat  towards  the  seasonal  mean  the  upper  levels  reach  a  relative 
minimum.  Finally,  the  Min-to-Max  inflection  point  (Trace  4)  completes 
the  cycle  as  the  upper-level  heights  shift  towards  the  seasonal  average 
and  the  lower-level  heights  become  positive.  This  lag  between  lower  and 
upper  levels  is  significant  because  it  implies  a  vertical  structure  exists  in 
the  oscillation  where  the  lower  levels  lead  the  upper  levels.  This  implied 
vertical  tilt  suggests  the  oscillation  could  have  a  baroclinic  structure. 

To  summarize,  it  appears  the  phase  of  the  12-d  oscillation  during  dry  season 
1984  can  be  detected  using  the  vertical  composites  of  all  three  variables.  Although 
the  composite  anomalies  may  be  too  small  to  determine  the  phase  of  the  oscillation 
using  quantitative  methods,  the  trends  of  the  anomalies  appear  to  be  consistent  and 
well-defined  throughout  the  cycle  of  the  oscillation,  especially  the  v-wind  component 
and  height  composites.  These  trends  would  allow  forecasters  to  monitor  and  forecast 
each  phase  of  the  oscillation. 
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Wet  Season  1984  Vertical  Composites 

The  composites  for  wet  season  were  constructed  using  the  same  methods  as 
the  dry  season.  The  results  of  the  IMSL  periodogram  analyses  (not  shown)  are 
summarized  as  follows: 

•  Very  broad  spectral  peaks  in  the  30-60  d  range  were  detected  at  all  levels 
for  the  u-wind  component.  Spectral  peaks  corresponding  to  the  period 
of  15-25  d  were  found  at  some  levels.  The  strongest  and  most  persistent 
spectral  peaks  were  in  the  11-12  d  range  and  were  found  at  all  levels,  but 
were  strongest  in  the  upper  levels. 

•  Strong  spectral  peaks  were  found  in  the  11-12  d  range  at  all  levels  for  the 
v-wind  component.  Some  concentration  of  spectral  power  was  detected 
at  several  levels  in  both  the  30-60  d  and  approximately  15-25  d  periods. 

•  Very  strong  spectral  peaks  in  the  height  fields  were  detected  at  around 
50-60  d  in  the  upper  levels  along  with  weak  spectral  peaks  in  the  15-30  d 
range.  Persistent  peaks  in  the  11-13  d  period  were  found  at  all  levels  but 
were  strongest  in  the  upper  levels. 

The  IMSL  periodogram  analysis  results  for  wet  season  were  similar  to  those  of 
dry  season.  The  12-d  oscillation  was  detected  at  all  levels  in  all  three  variables,  but 
was  again  strongest  in  the  v-wind  component.  Next,  the  12-d  filter  was  applied  to 
each  variable  at  every  level  and  composite  diagrams  were  constructed  in  the  same 
manner  as  for  dry  season.  The  vertical  composites  and  a  brief  discussion  of  their 


trends  are  as  follows: 
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•  The  u-wind  composites  are  shown  in  Figs.  39,  40,  41,  and  42.  Very  latge 
deviations  from  the  seasonal  mean  can  be  seen  in  the  upper  levels  (20- 
15  kPa)  as  the  composite  diagrams  follow  the  12-d  cycle.  Beginning  with 
the  Maximum  composite  (Fig.  42),  the  seasonal  mean  shows  easterly  winds 
of  about  4-6  kt  in  the  20-15  kPa  layer,  but  the  composite  winds  show  a 
westerly  wind  component  of  around  2-3  kt  during  this  phase.  Moving 
ahead  3  d  to  the  Max-to-Min  inflection  point  composite  (Fig.  43).  the 
composite  winds  are  nearly  equal  to  the  seasonal  average.  However,  3  d 
later  in  the  Minimum  composite  (Fig.  44),  the  composite  u-component 
winds  are  now  easterly  at  over  10  kt.  This  represents  a  deviation  of  12- 
13  kt  in  the  u-wind  component  over  the  6  d  period  from  the  Maximum 
composite  to  the  Minimum  composite,  as  well  as  a  switch  from  weak 
westerlies  to  strong  easterlies. 

•  The  vertical  composites  for  the  v-wind  component  are  shown  in  Figs.  43, 
44,  45,  and  46.  Strong  directional  shifts  in  the  upper-level  v-wind  com¬ 
ponent  can  also  be  seen  in  the  vertical  composites.  Beginning  with  the 
Maximum  composite  (Fig.  46),  the  composite  winds  are  roughly  equal  to 
the  seasonal  average  in  the  upper  levels.  Moving  to  the  Max-to-Min  in¬ 
flection  point  (Fig.  47),  the  composite  wind  at  20  kPa  is  now  southerly  at 
approximately  5  kt.  Three  days  later  at  the  Minimum  composite  (Fig.  48), 
the  composite  winds  are  again  nearly  equal  to  the  seasonal  average,  but 
on  the  Min-to-Max  inflection  point  composite  the  winds  are  now  from 
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Figure  39.  Vertical  composite  for  u-wind  component  corresponding  to  Maximum 
OLR  days,  wet  season  1984.  Solid  line  represents  seasonal  average, 
dashed  line  represents  deviation. 
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Figure  40.  Vertical  composite  for  u-wind  component  corresponding  to  Max-to-Min 
OLR  inflection  points,  wet  season  1984.  Solid  line  represents  seasonal 
average,  dashed  line  represents  deviation. 
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(EASTERLY)  KNOTS  (WESTERLY) 


Figure  41.  Vertical  composite  for  u-wind  component  corresponding  to  Minimum 
OLR  days,  wet  season  1984.  Solid  line  represents  seasonal  average, 
dashed  line  represents  deviation. 
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Figure  42.  Vertical  composite  for  u-wind  component  corresponding  to  Min-to-Max 
OLR  inflection  points,  wet  season  1984.  Solid  line  represents  seasonal 
average,  dashed  line  represents  deviation. 
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(NORTHERLY)  KNOTS  (SOUTHERLY) 


Figure  43.  Vertical  composite  for  v-wind  component  corresponding  to  Maximum 
OLR  days,  wet  season  1984.  Solid  line  represents  seasonal  average, 
dashed  line  represents  deviation. 
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Figure  44.  Vertical  composite  for  v-wind  component  corresponding  to  Max-to-Min 
OLR  inflection  points,  wet  season  1984.  Solid  line  represents  seasonal 
average,  dashed  line  represents  deviation. 
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Figure  45.  Vertical  composite  for  v-wind  component  corresponding  to  Minimum 
OLR  days,  wet  season  1984.  Solid  line  represents  seasonal  average, 
dashed  line  represents  deviation. 
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Figure  46.  Vertical  composite  for  v-wind  component  corresponding  to  Min-to-Max 
OLR  inflection  points,  wet  season  1984.  Solid  line  represents  seasonal 
average,  dashed  line  represents  deviation. 
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the  north  at  roughly  5-6  kt.  So  again,  within  the  span  of  6  d  the  com¬ 
posite  v-component  winds  have  changed  directions  from  a  5  kt  southerly 
component  to  a  5-6  kt  northerly  component. 

•  Figure  47  shows  the  vertical  composites  of  the  height  anomalies.  The 
height  deviations  can  again  be  detected  in  both  the  upper  and  lower 
levels,  similar  to  the  dry  season  height  composites.  However,  this  time  the 
upper- level  height  anomalies  lag  the  lower-level  height  anomalies,  as  can 
be  seen  from  following  the  cycle  of  the  12-d  oscillation  in  the  composites. 
This  structure  is  opposite  from  that  found  in  the  dry  season  and  would 
imply  that  a  different  mechanism  is  responsible  for  the  OLR  anomalies 
associated  with  a  period  of  12  d. 

The  composites  for  wet  season  show  that  strong  trends  in  all  three  variables  are 
associated  with  the  different  phases  of  the  12-d  oscillation, with  the  wet  season  trends 
are  noticeably  stronger  than  the  dry  season  trends.  The  next  section  will  attempt  to 
describe  the  synoptic  conditions  that  could  be  associated  with  the  12-day  oscillation 
using  the  information  obtained  from  the  vertical  composites. 

Possible  Synoptic  Scenarios  Associated  with  the  12-Day  Oscillation 

Dry  Season  1984 

A  summary  of  the  resultant  wind  field  and  height  anomalies  for  each  phase  of 
the  12-d  oscillation  during  dry  season  1984  is  shown  in  Table  1.  The  data  in  this 
table  were  obtained  from  the  vertical  composites  presented  in  the  previous  section. 


KPA 
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Figure  47.  Vertical  composite  for  height  anomalies  for  wet  season  1984.  Solid  line 
along  zero  deviation  line  represents  seasonal  mean,  sol’d  uaces  represent 
deviations  from  the  mean.  Trace  1  corresponds  to  Maximum  OLR 
composite,  Trace  2  to  Max- to- Min  composite.  Trace  3  to  Minimum  OLR 
composite,  Trace  4  to  Min-to-Max  composite. 
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Table  1.  Resultant  winds  and  height  anomalies  for  Panama  corresponding  to  the 
four  phases  of  the  12-d  oscillation  for  dry  season,  1984.  The  resultant  wind 
speeds  are  in  knots. 


25  I\PA 


85  KPA 


MAX 

MAX-MIN 

MIN 

MIN-MAX 

HEIGHT 

+ 

0 

— 

0 

WIND 

250/03 

210/07 

175/10 

170/06 

HEIGHT 

0 

— 

0 

+ 

WIND 

020/14 

015/11 

015/12 

015/14 

The  height  anomalies  are  depicted  as  relative  maxima  (4-)  and  minima  (-).  while 
zero  height  anomalies  indicate  the  height  composi+e  values  were  near  the  seasonal 
mean  for  that  phase.  The  selection  of  single  levels  to  represent  the  upper  and  lower 
troposphere  was  difficult  in  that  the  strongest  anomalies  for  each  variable  we'  ’  not 
always  found  at  the  25  kPa  and  85  kPa  levels.  For  example,  during  wet  season  the 
strongest  shift  in  the  v-wind  component  anomalies  is  found  at  the  20  kPa  level,  with 
a  weaker  signal  at  25  kPa.  However,  the  forecaster  trying  to  predict  the  phase  of  the 
12-d  oscillation  will  not  have  the  benefit  of  having  the  current  season  u-and  v-wind 
component  and  height  anomaly  composites  to  know  which  level  contains  the  strongest 
signal.  It  is  also  very  likely  that  the  structure  of  these  composites  will  vary  a  great 
deal  from  year-to-year.  With  that  in  mind,  the  25  and  85  kPa  levels  were  chosen 
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because  they  consistently  capture  the  trend,  and  usually  the  maximum  anomalies,  of 
the  deviations  found  throughout  the  12-d  cycle  for  both  dry  and  wet  seasons. 

Starting  the  12-d  cycle  shown  in  Table  1  with  the  Maximum  phase  (correspond¬ 
ing  to  maximum  OLR  values  over  Panama),  positive  height  anomalies  are  found  at 
25  kPa  while  the  composite  heights  are  near  the  seasonal  mean  at  85  kPa.  Moving 
through  the  four  phases  of  the  cycle  it  can  a ra;  e  een  that  the  lower  level  height 
anomalies  lead  those  of  the  upper-levels,  suggesting  a  vertical  tilt  to  the  oscillation. 
The  wind  composites  show  that  the  majority  of  the  fluctuation  in  the  wind  field  are 
found  in  the  upper- level  winds,  with  the  lower- level  wind  flow  remaining  fairly  con¬ 
stant  as  the  strong  northeasterly  winds  typical  of  the  dry  season  dominate.  Knowing 
this,  we  will  now  concentrate  on  the  fluctuations  contained  in  the  upper-level  winds. 
Fig.  48  shows  the  upper-level  wind  flow  from  climatology  for  the  month  of  January 
(Sadler  (1975)),  which  is  representative  of  dry  season  conditions.  The  main  features 
around  the  Central  American  region  are  a  weak  ridge  over  Panama  and  the  rest  of 
Central  America  with  a  broad  trough  located  to  the  west.  Relating  the  upper-level 
wind  composites  of  Table  1  to  climatology,  the  following  scenario  is  suggested  (see 
Fig.  49): 

•  During  the  Maximum  phase,  the  upper-level  trough  usually  located  to 
the  west  of  Panama  has  moved  to  the  east  of  the  country,  resulting  in 
clear  conditions  (positive  OLR  anomalies)  over  the  region.  Weak  west- 
southwesterly  flow  over  Panama  and  positive  height  anomalies  are  the 
result  of  a  broad  ridge  aloft  just  to  the  east  of  the  region. 


Figure  48.  20  kPa  streamline  climatology  for  January 


Maximum  2.  Max- 


e  49.  Possible  synoptic  scenario  associated  with  the  four  phases  of  the  12-d  oscillation,  dry  season 
represent  OLR.  anomalies  greater  than  i  (i  Win  2. 


86 


•  Three  days  later,  a  broad  trough  begins  to  build  far  to  the  west  of  Central 
America,  with  strong  negative  OLR  anomalies  developing  just  ahead  of  it. 
Southeasterly  winds  aloft  and  the  decrease  of  the  height  anomalies  over 
Panama  suggest  the  upper-level  ridge  has  moved  off  to  the  east. 

•  At  the  Minimum  phase  of  the  cycle  the  negatively-tilted  trough  has  now 
moved  just  off  the  west  coast  of  Central  America  and  has  deepened,  as 
indicated  by  the  relatively  strong  south-southeast  upper-level  flow  and 
the  negative  height  anomalies  over  Panama.  Note  that  the  strongest 
convective  activity  (shown  by  the  negative  OLR  anomalies)  is  also  found 
over  Panama  and  Central  America,  just  ahead  of  this  trough. 

•  Three  days  later  at  the  Min-to-Max  inflection  point  it  appears  the  trough 
has  weakened  and  moved  just  to  the  northeast  of  Panama,  as  reflected 
by  the  weakened  south-southeast  flow  and  an  increase  in  the  upper-level 
height  anomalies  from  a  relative  negative  to  the  seasonal  mean.  At  this 
point  the  cycle  repeats  again  as  weak  ridging  begins  to  build  over  the 
eastern  Pacific. 

While  it  is  unrealistic  to  expect  to  be  able  to  analyze  the  upper-level  wind  field 
of  the  entire  domain  using  winds  from  only  one  location,  this  scenario  does  agree  with 
climatology  and  the  location  of  the  OLR  anomalies,  and  also  provides  an  explanation 


for  the  apparent  west-to-east  movement  of  the  OLR  anomalies. 
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Table  2.  Resultant  winds  and  height  anomalies  for  Panama  corresponding  to  the 
four  phases  of  the  12-d  oscillation  for  wet  season,  1984.  The  resultant  wind 
speeds  are  in  knots. 


MAX 


25  KPA 


HEIGHT 


WIND 

120/05 

HEIGHT 

0 

WIND 

045/07 

MAX-MIN 

0 


140/07 


040/06 


MIN 

MIN-MAX 

+ 

0 

100/09 

080/08 

0 

+ 

025/04 

035/05 

Wet  Season  1984 

Table  2  shows  a  summary  of  the  resultant  winds  and  height  anomalies  for  wet 
season  1984.  The  data  were  obtained  in  the  same  manner  as  dry  season.  As  was  the 
case  in  the  dry  season,  the  majority  of  the  variation  in  the  wind  field  is  contained  in 
the  upper  levels  (25  kPa),  as  the  lower-level  winds  (85  kPa)  remain  fairly  constant 
from  the  northeast  throughout  the  12-d  cycle.  One  major  difference  between  the  two 
seasons  is  that  the  upper-level  height  anomalies  appear  to  lead  those  in  the  lower 
levels,  which  is  opposite  from  the  dry  season  pattern.  An  examination  of  the  upper- 
level  wind  climatology  for  August  shown  in  Fig.  50  reveals  that  the  controlling  feature 
for  Panama  appears  to  be  the  large  anticyclone  associated  with  the  subequatorial 
ridge  located  to  the  northeast  of  Panama.  The  resultant  wind  field  of  Table  2.  along 
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with  climatology,  suggests  that  the  following  synoptic  sequence  (shown  in  Fig.  51) 
could  be  associated  with  the  12-d  oscillation: 

•  Beginning  with  the  Maximum  phase,  the  anticyclone  that  is  shown  in 
climatology  is  located  well  to  the  northeast  of  Panama,  resulting  in  weak 
southeasterly  winds  and  a  negative  pressure  anomaly.  The  outflow  aloft 
associated  with  this  anticyclone  is  also  well  to  the  north  of  the  region, 
resulting  in  relative  maximum  OLR  values  over  Central  America. 

•  Three  days  later  at  the  Max-to-Min  inflection  point,  the  anticyclone  begins 
to  move  slowly  to  the  west.  This  is  reflected  by  the  increase  in  convective 
activity  over  Panama  and  the  continued  southeasterly  flow. 

•  During  the  Minimum  phase  of  the  12-d  oscillation  the  anticyclone  has 
continued  to  move  to  the  west  and  appears  to  be  located  just  to  the  north 
of  Panama,  as  shown  by  the  upper-level  winds  which  are  now  from  the 
east-southeast.  At  this  time  the  maximum  outflow  aloft  associated  with 
the  anticyclone  and  the  subequatorial  ridge  is  located  over  Panama  and 
Central  America,  resulting  in  the  large  area  of  negative  OLR  anomalies 
over  the  region  which  extends  along  the  monsoon  trough  to  the  west. 

•  By  the  final  phase  of  the  oscillation,  the  Min-to-Max  inflection  point, 
the  anticyclone  appears  to  be  drifting  slowly  to  the  west-northwest  as 
Panama’s  upper-level  winds  are  now  from  the  east-northeast.  The  mini¬ 
mum  OLR  anomalies  associated  with  the  anticyclone  have  also  drifted  to 
the  west-northwest  and  are  located  off  the  coasts  of  Guatemala  and  El 
Salvador.  At  this  point  it  is  not  clear  from  the  information  available,  but 
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it  appears  the  anticyclone  dissipates  once  it  moves  to  the  west  and  another 
forms  to  the  northeast  of  Panama  and  the  cycle  begins  again.  Again,  these 
scenarios  are  difficult  to  verify  using  data  from  only  one  location. 

To  summarize,  the  synoptic  scenario  for  the  wet  season  depicted  above  differs 
from  that  of  dry  season.  In  this  case  the  factor  controlling  the  development  and 
movement  of  the  OLR  anomalies  through  the  12-d  oscillation  appears  to  be  the 
westward  movement  of  the  anticyclone  which  is  shown  by  climatology  to  be  located  to 
the  northeast  of  Panama.  Also,  the  OLR  anomalies  appear  to  develop  over  Panama 
and  drift  to  the  west-northwest  during  the  wet  season  as  opposed  to  the  west-to-east 
movement  observed  during  dry  season. 

Case  Study 

The  previous  sections  have  concentrated  on  using  upper-air  data  from  a  single 
station  to  detect  and  monitor  the  phase  of  the  12-d  oscillation.  This  section  will 
follow  a  complete  cycle  of  the  12-d  oscillation  during  wet  season  1984  through  the 
use  of  satellite  imagery.  The  satellite  imagery  is  from  the  Geostationary  Operational 
Environmental  Satellite  (GOES)  system.  The  satellite  imagery  shown  in  Figs.  52,  53, 
54,  and  55  are  enhanced  infrared  photos  with  7  km  resolution.  Clark  (1983)  describes 
the  color  enhancement  levels  as  follows: 

•  The  first  level  contours  (medium  and  light  gray)  represent  cloud  tops 
between  approximately  9,500  and  11,000  m. 

•  The  next  dark  gray  contour  indicates  thunderstorm  activity  with  maxi¬ 
mum  cloud  tops  of  approximately  12,500  m. 
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Figure  52.  GOES  enhanced  IR  satellite  photo  for  1800  GMT,  3  October  1984. 
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Figure  53.  GOES  enhanced  IR  satellite  photo  for  1900  GMT,  6  October  1984. 
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•  The  black  enhancement  represents  cloud  tops  to  approximately  13,400  m. 

•  The  gray  shade  embedded  within  the  black  contour  indicates  cloud  tops 
of  approximately  14,700  m. 

•  The  white  enhancement  within  the  black  contour  represents  overshooting 
cloud  tops  above  14,700  m. 

Using  this  enhancement,  the  thunderstorms  around  the  Panama  area  are  easily 
seen  in  Fig.  52.  This  photo  is  from  the  afternoon  of  3  October  1984,  in  the  middle  of 
the  wet  season.  This  day  was  shown  by  the  12-d  filtered  OLR  time  series  (Fig.  19)  to 
be  a  Minimum  OLR  day,  that  is,  a  day  with  more  convection  than  usual  in  the  Panama 
area.  Note  the  similarities  between  this  photo  and  the  Minimum  OLR  composite 
(Fig.  28).  The  satellite  photo  shows  strong  thunderstorm  activity  throughout  the 
Panama  region,  extending  to  the  west  along  the  monsoon  trough.  The  minimum 
OLR  composite  also  depicts  these  conditions  during  this  phase  of  the  12-d  oscillation. 
Note  the  similarities  at  around  15°N,  100°W  where  the  maximum  OLR  anomalies  of 
the  OLR  composite  correspond  with  the  relatively  clear  area  on  the  satellite  photo. 

Moving  ahead  3  d  to  Fig.  53  (Min-to-Max  phase),  the  thunderstorm  activity 
over  the  Panama  area  has  decreased  in  intensity,  and  a  large  area  of  clouds  and 
thunderstorm  activity  is  centered  at  approximately  10°N,  95°W.  Comparing  this 
photo  with  the  Min-to-Max  composite  of  Fig.  29,  the  areas  of  activity  match  up 
almost  exactly. 

The  satellite  photo  depicting  the  next  phase  of  the  cycle  is  shown  in  Fig.  54, 
which  is  the  Maximum  OLR  phase.  At  this  time  the  Panama  region  appears  to  have 
a  fair  amount  of  high  cloud  cover  with  some  embedded  showers,  but  the  cloud  cover 
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and  thunderstorm  intensity  has  decreased  significantly  as  compared  to  the  Minimum 
and  Min-to-Max  phases.  Also,  the  rest  of  Central  America  is  not  experiencing  any 
significant  thunderstorm  activity  in  this  afternoon  photo  which  is  unusual  for  this 
time  of  year.  The  OLR  composite  of  Fig.  26  corresponds  to  this  phase  of  the  cycle, 
and  the  relatively  clear  conditions  over  Central  America  seen  on  the  satellite  photo 
are  depicted  on  the  composite.  One  area  of  thunderstorm  activity  on  the  satellite 
photo  is  centered  at  approximately  10°N,  100°W,  and  this  area  matches  well  with 
the  area  of  negative  OLR  anomalies  shown  on  the  composite.  Also  worth  noting  is 
the  lack  of  convection  along  the  monsoon  trough  as  compared  to  the  Minimum  OLR 
phase. 

The  final  phase  of  the  12-d  oscillation  is  represented  by  the  satellite  photo  in 
Fig.  55.  At  this  time  strong  thunderstorm  activity  is  beginning  to  develop  over  the 
Panama  area,  but  the  rest  of  Central  America  remains  relatively  cloud-free.  These 
general  conditions  also  show  up  on  the  composite  diagram  of  Fig.  27. 

The  variations  in  convection  and  cloud  cover  associated  with  this  12-d  oscillation 
may  not  always  be  as  clear-cut  as  this  case  study.  However,  by  following  the  trends  of 
the  upper-air  analyses  outlined  in  this  chapter  and  maintaining  continuity  of  satellite 
imagery,  the  cycle  of  the  12-d  oscillation  car,  be  followed.  This  would  allow  for 
forecasts  beyond  24  h  to  be  made  for  above/below  normal  amounts  of  convection  and 
cloud  cover. 

Correlation  Between  OLR  Data  and  Precipitation 

This  study  has  concentrated  on  determining  the  variability  of  convective  activity 
over  the  Panama  region  through  use  of  OLR  data.  One  additional  question  that 
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remains  concerns  the  ability  of  OLR  data  to  infer  actual  rainfall  totals.  If  low  OLR 
values  are  related  to  convective  activity,  then  there  should  be  a  high  correlation 
between  OLR  and  rainfall  amounts.  The  long-period  seasonal  trends  were  removed 
from  the  time  series  of  daily  precipitation  totals  from  raingauges  in  the  Panama 
Canal  region  (see  Fig.  3).  This  procedure  was  similar  to  that  used  on  the  daily 
average  OLR  time  series.  Once  rainfall  data  were  in  this  form,  the  degree  of  linear 
association  between  the  rainfall  data  and  the  OLR  time  series  of  Figs.  9,  10,  11, 
and  12  was  calculated  using  the  coefficient  of  determination  (Neter,  Wasserman,  and 
Kutner  (1983)),  where 


2  =  (  E (X,ZXM-Y)_ 

The  results  of  these  calculations  are  as  follows: 

•  Dry  season  1984:  r2=.05 

•  Wet  season  1984:  r2=.07 

•  Dry  season  1985:  r2=.01 

•  Wet  season  1985:  r2=.13. 


(11) 


The  calculations  were  again  performed,  this  time  using  only  the  two  OLR  points 
from  the  domain  (shown  in  Fig.  6)  that  were  closest  to  Panama.  Again,  these  two 
OLR  points  were  averaged  into  one  daily  OLR  value  for  the  area  and  the  seasonal 
trends  were  removed.  The  coefficient  of  determination  values  relating  this  smaller 


area  and  Panama  rainfall  are  as  follows: 
•  Dry  season  1984:  r2=.08 
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•  Wet  season  1984:  r2=.10 

•  Dry  season  1985:  r2=.04 

•  Wet  season  1985:  r2=.13 

The  values  of  r2  can  range  from  0  (indicates  no  linear  association  between  the 
two  data  sets)  to  1  (indicates  strong  linear  association).  The  values  obtained  here 
would  seem  to  imply  that  there  is  little  correlation  between  the  OLR  averages  of 
either  area  and  Panama  rainfall.  These  numbers  are  significantly  lower  than  what 
was  expected,  and  it  was  thought  that  one  possible  cause  for  the  low  r2  values  could 
be  traced  to  the  differences  in  the  techniques  used  to  measure  the  two  variables.  The 
rainfall  data  are  comprised  of  24  hourly  totals  from  the  26  raingauge  sites  which 
record  rainfall  continually  throughout  the  observation  period.  The  OLR  data  are 
different  in  that  they  only  represent  snapshots  of  conditions  taken  twice  a  day,  at 
approximately  0230  and  1430  LST.  So  it  is  doubtful  that  the  OLR  data,  having  only 
those  two  observation  times,  would  coincide  with  the  daily  rainfall  periods.  One 
solution  was  to  use  only  the  0200-0300  and  1400-1500  LST  rainfall  totals  to  correlate 
with  the  OLR  data.  To  eliminate  any  bias  concerning  the  spatial  extent  of  the  rain 
resulting  from  one  raingauge  recording  a  downpour  of  several  inches  in  an  hour  while 
the  other  raingauges  were  dry,  the  frequency  of  precipitation  was  recorded.  If  any 
rain  fell  at  a  raingauge  site  during  the  specified  hour,  regardless  of  intensity,  a  T’  was 
recorded,  while  if  no  rain  was  measured  a  ‘0’  was  recorded.  For  example,  if  6  of  the  26 
total  raingauge  sites  recorded  precipitation  during  0200-0300  LST  and  14  of  the  26 
raingauges  recorded  precipitation  during  the  1400-1500  LST  period,  the  daily  rainfall 
total  used  to  correlate  with  the  OLR  daily  average  would  be  20  (out  of  a  total  of  52 
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raingauge  measurements).  The  coefficient  of  determination  was  again  calculated  for 
the  both  OLR  domains  (the  20-point  area  average  as  well  as  the  2-point  average), 
however,  the  results  (not  shown)  were  not  significantly  higher. 

There  could  be  many  factors  contributing  to  the  low  r2  values,  some  of  which 
could  be: 

•  The  OLR  averages  used  in  this  study  are  measured  and  averaged  over  a 
much  larger  spatial  scale  than  that  of  the  precipitation.  Even  the  smaller 
area  of  OLR  used  in  the  correlation  study  was  the  average  of  two  2.5° 
latitude-longitude  areas.  Local  conditions,  such  as  orographically-induced 
precipitation  that  could  be  causing  isolated  precipitation  over  the  Panama 
Canal  region,  would  not  be  reflected  on  this  large  scale. 

•  Some  of  the  precipitation,  especially  over  the  highlands  of  central  Panama, 
could  be  occurring  with  low,  warm  cloud  tops  having  relatively  high 
OLR  values.  This  would  contradict  the  original  assumption  that  all 
precipitation  over  the  Panama  area  is  associated  with  cold  cloud  tops 
associated  with  intense  convection. 

•  Both  the  OLR  and  precipitation  time  series  are  extremely  noisy.  A 
smoothing  function  applied  to  the  time  series  may  reduce  some  of  the 
erratic  fluctuations  that  could  be  contributing  to  a  low  r2  value. 

In  conclusion,  it  appears  in  this  study  that  the  OLR  data  set  can  not  be  used 
to  directly  predict  rainfall  events  in  the  Panama  Canal  area.  Instead,  the  OLR  data 
depict  fluctuations  in  cloud  cover  over  the  region,  and  these  variations  imply  greater 
than  or  less  than  normal  chances  of  precipitation  over  the  region. 
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CHAPTER  VII 

MODULATION  OF  12-DAY  OSCILLATION 

In  Chapter  VI  it  was  mentioned  that  the  12-d  filtered  OLR  time  series  for  wet 
season  1984  (Fig.  19)  appeared  to  contain  a  modulation  of  the  signal  on  the  time 
scale  of  70-80  d.  While  looking  through  daily  satellite  imagery  for  case  study  photos 
an  interesting  relationship  between  the  “active”  and  “inactive”  phases  of  the  12-d 
oscillation  and  tropical  storm  activity  in  eastern  Pacific  was  discovered.  Figure  56  is 
the  12-d  filtered  OLR  time  series  for  wet  season  1984  (which  is  the  same  as  Fig.  19), 
which  has  been  subjectively  divided  into  active  and  inactive  phases.  Active  phases 
are  periods  when  the  12-d  signal  is  strongest,  while  inactive  phases  contain  relatively 
weak  12-d  signals.  The  black  bars  running  along  the  horizontal  axis  show  the  duration 
in  days  of  tropical  systems  that  developed  in  the  eastern  Pacific  (Gunther  and  Cross 
(1985)).  Note  the  large  number  of  systems  that  existed  during  the  inactive  phases 
as  compared  to  the  active  phases.  Quantitatively,  the  total  number  of  systems  that 
developed  during  the  inactive  phases  was  16  (5  tropical  depressions,  2  tropical  storms, 
and  9  hurricanes)  compared  to  only  8  (3  tropical  depressions,  3  tropical  storms,  and 
2  hurricanes)  for  the  active  phase.  A  slight  modulation  of  the  1985  wet  season  12-d 
filtered  time  series  (Fig.  21)  can  be  seen,  but  it  is  not  nearly  as  strong  as  that  of  1984 
wet  season. 

The  reasons  for  this  connection  between  the  relative  strength  of  the  12-d 
oscillation  and  the  development  of  tropical  storms  in  the  eastern  Pacific  are  not 
apparent  from  this  research.  Possibly  a  north-south  oscillation  of  the  subequatorial 
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lines  along  horizontal  axis  represent  tropical  cyclone  activity. 
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ridge  exists  during  the  wet  season  which  creates  favorable  and  unfavorable  upper- 
level  environments  for  tropical  storm  development  in  the  eastern  Pacific.  During  the 
inactive  phase  when  favorable  conditions  exist  in  the  eastern  Pacific  for  tropical  storm 
development,  the  upper-level  anticyclone  associated  with  the  subequatorial  ridge  (see 
Figs.  50  and  51)  could  be  too  far  to  the  north  of  Panama.  This  could  result  in  a  weaker 
12-d  oscillation  by  removing  the  outflow  mechanism  thought  to  be  responsible  for  the 
periodic  fluctuations  in  convection.  However,  significantly  more  research  needs  to  be 
done  on  this  subject  before  any  theories  can  be  proven. 
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CHAPTER  VIII 

SUMMARY  AND  AREAS  FOR  FURTHER  STUDY 

Convective  variability  over  the  Panama  region  was  examined  through  spectral 
analysis  of  OLR  time  series  data.  The  strongest  and  most  persistent  spectral  peak 
was  found  in  the  12-d  range,  with  weaker  spectral  power  evident  in  the  30-60  d  and 
15-30  d  ranges.  The  OLR  time  series  for  the  Panama  region  was  filtered  for  12  d,  and 
composite  diagrams  corresponding  to  each  of  the  four  phases  (Maximum,  Minimum, 
and  two  inflection  points)  of  the  12-d  oscillation  were  constructed  for  the  domain. 
The  following  conclusions  were  drawn  from  the  composite  analyses: 

•  The  composites  for  dry  season  1984  revealed  west-to-east  propagating 
OLR  anomalies  that  appear  to  originate  in  the  eastern  Pacific  and  build 
to  the  east  over  Panama. 

•  The  composites  for  wet  season  1984  differed  significantly  from  dry  season 
in  that  the  OLR  anomalies  appear  to  originate  over  the  Panama  region 
and  drift  slowly  to  the  west -northwest. 

Vertical  composites  of  upper-air  data  from  Albrook  Air  Force  Station  were 
then  constructed  to  correspond  to  each  of  the  four  phases  of  the  12-d  oscillation. 
Composites  of  u-  and  v-wind  components  and  height  of  the  pressure  surfaces  were 
examined  in  order  to  determine  the  vertical  structure  of  the  oscillation.  The  vertical 


composites  revealed  the  following: 
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•  In  dry  season  1984,  12-d  fluctuations  from  the  seasonal  mean  were  evident 
in  both  the  u-  and  v-wind  components  in  the  upper-levels,  with  the  v- 
wind  component  showing  the  strongest  fluctuations.  Cyclic  fluctuations 
were  present  in  the  height  composites  at  all  levels  during  the  dry  season. 
These  height  fluctuations  reflected  a  baroclinic  structure  to  the  oscillation 
as  the  lower-level  height  anomalies  led  the  upper-level  height  anomalies 
by  approximately  a  quarter-cycle. 

•  Fluctuations  of  12-d  were  again  evident  in  the  upper-level  u-  and  v-wind 
composites  and  at  all  levels  of  the  height  composites  during  the  wet  season. 
These  fluctuations  were  significantly  stronger  than  those  detected  in  the 
dry  season.  The  major  difference  between  the  two  seasons  was  found  in 
the  height  composites.  These  composites  suggest  the  vertical  structure 
of  the  oscillation  is  reversed  in  the  wet  season  as  the  upper-level  height 
anomalies  now  lead  the  lower-level  height  anomalies  by  approximately  a 
quarter-cycle. 

The  fluctuations  of  the  departures  of  the  three  variables  from  their  seasonal 
means  would  appear  to  be  difficult  to  track  from  a  purely  quantitative  standpoint, 
especially  since  the  seasonal  mean  for  the  season  in  progress  can  not  be  determined. 
However,  the  trends  of  each  variable  through  each  of  the  four  phases  appear  to  be 
strong  enough  to  enable  forecasters  to  monitor  the  progress  of  the  12-d  oscillation 
and  the  associated  variations  in  convective  activity. 
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Using  OLR  and  upper-air  composites  as  well  as  upper-level  wind  climatology, 
synoptic  scenarios  associated  with  the  12-d  oscillation  were  constructed  for  each 
season.  The  results  are  as  follows: 

•  The  upper-level  climatological  winds  show  predominantly  westerly  flow 
with  a  broad  trough  embedded  in  this  flow,  located  to  the  west  of  Central 
America  in  the  eastern  Pacific.  The  vertical  composites  suggest  that  this 
trough  periodically  deepens  and  moves  to  the  east  of  Panama,  resulting 
in  enhanced  convective  activity  ahead  of  it. 

•  The  wet  season  pattern  is  significantly  different  from  that  of  dry  season. 
During  this  time,  climatology  shows  the  upper-levels  over  Panama  to  be 
dominated  by  the  subequatorial  ridge  and  an  associated  anticyclone  to  the 
northeast  of  Panama.  The  composites  show  indications  that  this  anticy¬ 
clone  moves  slowly  to  the  west -northwest,  and  the  increased  outflow  aloft 
associated  with  this  anticyclone  enhances  convective  activity  over  Panama 
as  it  passes  just  north  of  the  region.  This  convection  associated  with  the 
anticyclone  continues  appears  to  drift  slowly  to  the  west-northwest  until 
the  anticyclone  weakens,  then  it  seems  another  anticyclone  develops  to 
the  northeast  of  Panama  and  the  cycle  begins  again. 

Another  interesting  feature  associated  with  this  12-d  oscillation  is  the  apparent 
link  between  tropical  cyclone  activity  in  the  eastern  Pacific  and  70-80  d  fluctuations 
in  the  strength  of  the  12-d  oscillation  over  Panama.  In  the  wet  season,  during  the 
periods  when  the  12-d  oscillation  over  Panama  was  strong  (active  periods),  tropical 
cyclone  activity  over  the  eastern  Pacific  was  relatively  weak  (3  tropical  depressions, 
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3  tropical  storms,  and  2  hurricanes).  By  comparison,  during  the  periods  when  the 
12-d  oscillation  over  Panama  was  weak  (inactive  periods),  a  significant  increase  in 
eastern  Pacific  tropical  cyclone  activity  (5  tropical  depressions,  2  tropical  storms,  and 
9  hurricanes)  was  noted.  Although  the  cause  for  this  modulation  was  not  addressed 
in  this  research,  one  possible  cause  could  be  periodic  north-south  oscillations  in  the 
position  of  the  subequatorial  ridge.  A  northward  migration  of  the  ridge  could  result 
in  a  weaker  12-d  signal  by  weakening  the  outflow  mechanism  necessary  for  strong 
convection  over  Panama,  while  at  the  same  time  creating  favorable  conditions  for 
tropical  cyclone  development  further  north  over  the  eastern  Pacific  where  the  low- 
level  monsoon  trough  lies. 

This  research  discussed  the  12-d  oscillation  and  its  effects  on  convective  activity 
over  the  Panama  region.  There  are  still  many  aspects  of  the  12-d  oscillation  to  be 
examined.  They  are: 

•  Examine  several  years  of  OLR  data  to  determine  the  interannual  variabil¬ 
ity  of  the  oscillation.  The  12-d  signal  was  very  strong  in  1984,  but  was 
not  quite  as  strong  in  1985.  Perhaps  some  general  circulation  feature  (El 
Nino,  quasi-biennial  oscillation)  is  connected  with  it. 

•  Although  not  shown,  the  1985  OLR  composites  showed  fairly  strong 
connections  between  Panama  and  portions  of  South  America  with  a  period 
of  12  d.  If  several  more  years  of  data  were  examined  it  could  be  determined 
if  this  was  just  a  one-time  occurrence  or  if  a  connection  does  exist  between 
the  northern  and  southern  hemispheres. 
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•  The  dry  season  OLR  composites  showed  the  OLR  anomalies  appeared 
to  extend  to  the  west  from  beyond  the  western  edge  of  the  domain  at 
120°W.  The  domain  could  be  expanded  to  the  west  to  determine  the 
western  extent  of  the  oscillation. 

•  A  more  detailed  study  of  the  upper-level  features  associated  with  this 
oscillation  could  be  performed  by  using  hemispheric  wind  data.  As  was 
noted,  it  is  difficult  to  imply  conditions  over  the  entire  domain  from 
only  one  vertical  sounding.  A  hemispheric  wind  study  could  confirm  the 
synoptic  patterns  inferred  in  this  research. 

•  A  closer  examination  could  be  made  of  the  relationship  between  tropical 
cyclone  activity  in  the  eastern  Pacific  and  the  relative  strength  of  the  12-d 


oscillation  over  Panama. 
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